Fermi-surface pockets in YBa2Cu306.5 : A comparison of ab initio techniques 
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We study the Fermi surface of metallic, non-magnetic ortho-11 YBa2Cu306.5 using three different 
density-functional-based band-structure techniques (GGA, GGA-I-U, PSIC). The calculated Fermi 
surface exhibits no pockets in GGA-I-U and PSIC, a minor one in GGA. Upon shifting the Fermi level 
in the vicinity of the calculated value, we instead observe several pocket structures. We calculate 
their cross-sectional areas and cyclotron masses. Overall, our calculations show no solid evidence of 
the existence of electron-like -nor, in fact, of any- Fermi surface pockets in this phase. This suggests 
that the origin of the pockets should be sought for in other, different phases. 

PACS numbers: 71.18.-|-y,74.72.-h,74.25. Jb 



I. INTRODUCTION 

The Fermi surface of underdoped high-temperature 
cuprate superconductors is currently under intense inves- 
tigation. Recentlj^i2ii^i^ Shubnikov-de Haas (SdH) and 
de Haas-van Alphen (dHvA) oscillations were observed 
in ortho-ll YBaaCugOe.s (henceforth YBCO; ortho-U 
stands for the chain-ahgned oxygen configuration with 
one Cu(l)-0 chain per 2x1x1 cell). These oscillations 
(of resistance and HaU coefficient in SdH, and magneti- 
zation in dHvA) correspond to closed sections (pockets) 
of the Fermi surface and they exhibit, as a function of 
the inverse of the magnetic field, characteristic frequen- 
cies related to the cross-sectional area of the pocket (or 
pockets: their number and location is undetermined). 

The frequency measured by dHvA experiments (more 
accurate than SdH) is 540±4 T, corresponding to a small 
portion (2%) of the Brillouin zone being enclosed by the 
pockets. The cyclotron mass, deduced from a Lifshitz- 
Koshevic fit of the oscillation amplitude vs temperature, 
is m=1.76±0.07 free-electron masses. The oscillations 
were observedi*^*^ in high field and only at low (4 K) tem- 
perature; the sign of the Hall coefficient was seen to be- 
come negative from about 25 K downward, and this was 
interpreted as a signature of the pockets in question being 
electron-like in nature. A further recent measurement^ 
in YBa2Cu306.5i reported, in addition to the same sig- 
nal of Refill, an oscillation with frequency and mass in 
the vicinity of 1600 T and 3.4 mg respectively, allegedly 
(see Ref|3, p. 201) associated with a hole-like pocket. 

SdH and dHvA examine at low temperature a state ob- 
tained by applying a high magnetic field to the supercon- 
ductor. To a first approximation this state is supposed to 
be the normal (possibly pseudo-gap) state. The simplest 
hypothesis is that once superconductivity is removed, 
YBCO is a metallic and non-magnetic Fermi-liquid like 
system (although more sophisticated options also exist, 
such as e.g. magnetic fluctuations and polaronic forma- 
tions in stripe-like morphology^ and more). Since exper- 
iments are often interpreted based on this assumption. 



an issue to be settled is whether or not the Fermi sur- 
face of this specific non-magnetic metallic phase exhibits 
pockets as revealed in experiments. If no calculated pock- 
ets exist, or can be identified with those observed, then 
some other phase will have to be invoked as the state 
accessed in oscillation experiments. To address this is- 
sue, here we employ three distinct techniques based on 
density- functional theory (DFT) : GGA (generalized gra- 
dient approximation), GGA-I-U, and PSIC (pseudo-Self- 
Interaction-Correction) method. Furthermore, we adopt 
the common practice (discussed below) of applying rigid- 
band shifts to explore the Fermi surface in a wide energy 
interval surrounding the calculated Ep. Our calculation 
widen the scope of recent^ calculations limited to the 
GGA approach. 

Our study shows that overall there is no reliable in- 
dication that non-magnetic metallic YBCO possesses 
electron-Fermi surface pockets. Specifically, only one 
technique (the GGA+U) finds an electron-like pocket, 
appearing however at a -60 meV shift away from calcu- 
lated Ep. None of the other techniques find any such 
pocket in a ±100 meV interval around Ep. As we will 
argue, in fact, there is only scant evidence for hole-like 
pockets as well. 

While we do not question the reliability of the 
SdH/dHvA experiments, we note that one may envis- 
age ways to generate a negative Hall coefficient other 
than the existence of electron-like pockets. For exam- 
ple, oscillations may be due to hole-like pockets, and the 
Hall coefficient positive-to-negative crossover may stem 
from the contribution of other electron-like Fermi sur- 
face structures, mixed up by differently temperature- 
dependent electron and hole mobilities. Other consider- 
ations that must be marked on the theoretical roadmap 
are that the pocket structure is partially at odds with 
the "Fermi arcs" observed^'^ '^"!^^'^^ in angle-resolved 
photo-emission spectroscopy; and that ordering phenom- 
ena, possibly related to magnetic structure or density 
waves, may be causing a reconstruction of the Fermi sur- 
face. 
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II. METHOD 

We calculate the band structure of YBCO in the 
non-magnetic metallic state with three different DFT- 
based techniques. We assumed the crystal struc- 
ture of YBa2Cu306.5 determined by Grybos et alM>^ 
We use the GGA (generalized gradient approxima- 
tion), GGA-I-U, and the pseudo-self-interaction correc- 
tion method (PSIC), a parameter- free, first-principles 
DFT-based methodic which correctly describes the 
physics of several correlated cuprateS )^^'^^i^^ and yet is 
practically viable for large-sized systems. In particular 
the PSIC is able to describe the competition of metallic 
and insulating phases of YBa2Cu306-i-a; from x—0 (where 
it iE^*> an antiferromagnetic Mott insulator) across two 
metal-insulator transitions^^ to metallic a;=l, obviating 
to the failures of plain GGA or similar approaches in this 
context. 

Our GGA and GGA-I-U calculations are carried out us- 
ing the VASP packagei2i2£ with the projector-augmented 
wave method (PAW)i2i The PSIC calculation are per- 
formed using a custom in-house code with ultrasoft 
pseudopotentials^^ and a plane wave basis set. The cutoff 
energy was set at 420 eV. A Monkhorst-Pack^^ 9x19x6 
grid was used for the self-consistency cycle. We inten- 
tionally used the in-plane 2x1 periodicity appropriate 
to chain-ordered ortho-ll YBCO at this specific dop- 
ing, since experiments are claimed to be performed in 
this structure. We tested non-spin-polarized calculations, 
spin-polarized calculations with small initial moments, 
and fixed-magnetic-moment calculations with zero im- 
posed magnetization, consistently getting the same re- 
sults, i.e. a non-magnetic metallic state. The Fermi sur- 
faces are visualized with the Xcrysden package. 

We used the Dudarev implementation^^ of GGA-I-U, 
whereby the relevant parameter is the difference U-J of 
the effective on-site Coulomb and exchange interactions. 
U-J was set to 9 eV for the d states of planar Cu, the 
value reproduces the fundamental gap of Mott-insulating 
antiferromagnetic YBa2Cu306 as obtained in PSI C^^i^^ 
or in experiment (no qualitative changes are observed 
down to f/-J=6 eV for YBCO). We underhne that the 
paramagnetic Fermi surface calculation is sensitive to [/- 
J via small orbital polarizations (i.e. deviations from ex- 
act half- filling) in the partially occupied Cu d^2_y2 states, 
and this may affect the details of band morphology in the 
vicinity of Ep. 



III. RESULTS AND DISCUSSION 



plays the bands within the GGA approach. Moving along 
the (7r/2,0)-(7r/2,7r) direction, the first band to cross Ep 
is mainly due to states of the Cu(l)-0 chain. This band 
is very close to being one dimensional. The next four 
bands crossings Ep come from the Cu02 planes. There 
is a splitting between the bonding and antibonding Cu02 
bands of -0.2 eV along the (7r/2,0)-(7r/2,7r) line at Ep. 
Each of these two bands are further split up by the ad- 
ditional 2a periodicity (this is most evident close to the 
point (7r/2,7r)). We find that the splitting is 40 meV at 
Ep along the (7r/2,0)-(7r/2,7r) direction. In the GGA cal- 
culation a fairly flat Cu(l)0 chain-Oapicai band crosses 
Ep close to the (0,7r) point and gives rise to a small tubu- 
lar quasi-2D hole pocket. This band is 13 meV above Ep 
at (0,7r). In addition, a second band with a similar char- 
acter lies just 20 meV below Ep at (0,7r). Our results are 
similar to calculations on YBCO reported previously.— 
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A. Band structures 
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FIG. 1: Band structures of YBa2Cu306.5. (a) GGA, (b) 
GGA+U with U-J^ 9 eV, (c) PSIC. 

In the central panel, Figlljb), we show the GGA+U 
bands. Overall, the GGA-I-U rendition appears quite 
close to those of GGA. This is expectable as U only af- 
fect magnetic and/or orbital-polarized states, thus the 
paramegnetic configuration is mildly affected. The main 
difference with respect to the GGA case is that the flat 
chain-apical bands crossing Ep near (0,7r) are now about 
80 meV above Ep and 140 meV below Ep at (0,7r), i.e. 
they are split by more than 200 meV, compared to about 
30 meV in GGA. This difference is due to the indirect 
(i.e. self-consistent) effect of the orbital polarization of 
in-plane Cu d states on the band manifoldi^ 

The right panel, FiglTJc), shows our calculation with 
the PSIC technique. Here we see more radical differences 
with respect to the other two methods, mainly due to the 
fact that PSIC corrects for sef-interaction Cu d as well 
as O p state occupations, so that the corrections can be 
equally sizable for non magnetic and/or non orbitally- 
polarized states. This description results in generally less 
dispersed band structure; chain bands are now far from 
Ep, and the net result is that there are no small pockets 
in the Fermi surface. 



In Fig[T]we compare the band structures in the k^^O 
plane, as obtained by the three methods. The dispersion 
in kz is weak and not important in the present context. 
kx and ky are in units of the inverse 1/a and 1/b of the in- 
plane lattice constants. The leftmost panel (Figl^a)) dis- 



B. Fermi surfaces 

Strictly speaking, the theoretical prediction of the 
Fermi surface is based on the calculated electronic struc- 



ture and Fermi level. Here, however, we also consider 
how the Fermi surface changes upon an upward or down- 
ward shift of the Fermi level compared to the calculated 
value. This is a fairly common practice in band theory 
studies of superconductors. The first motivation is that, 
while DFT calculations usually describe well the general 
features of the band structure of metals, small discrep- 
ancies in the relative positions of the bands are common 
when comparison with experiment is involved. (Gener- 
ally, this relates to structural details and of course to the 
DFT description of the electron correlation.) For exam- 
ple, in Sr2Ru04, studied in detail with the dHvA tech- 
nique, the Fermi energy needs to be shifted by 40 meV in 
either directional to improve the calculated-bands agree- 
ment with experiment. Even in MgB2, shifts of the order 
of 100 meV are needed.^^'-^ 

A further motivation pertaining to doped cuprates is 
that Fermi level shifts roughly simulate doping fluctu- 
ations. Of course the shift-doping relation depends on 
which specific band or bands are or get occupied upon 
shifting. In our case the maximum shifts applied (■~50-60 
meV) correspond to rather substantial doping fluctuation 
('^±0.04, i.e. a 30% of the nominal doping). 

In Fig[2]we collect the Fermi surface for the three tech- 
niques (top to bottom) , and upward to downward (left to 
right) shifts of the Fermi level. The top panel (Figl^Ja)) 
reports GGA results. For AEf=-I-50 meV the Fermi sur- 
face consists of just two large hole-like Cu02 sheets cen- 
tered on (7r/2,7r), plus three quasi-one-dimensional sheets 
(one from the chains, and two from the planes). As Ep 
shifts down, a small hole-like pocket develops near the 
(0,7r) point, originating from the flat CuO-Oapicai band 
discussed earlier. A further lowering of Ep causes this 
pocket to grow in size and then merge with the anti- 
bonding Cu02 plane sheet. As Ep is further reduced, 
the second CuO-BaO band crosses the Fermi level, giv- 
ing rise to another pocket. Eventually, this merges with 
the bonding Cu02 plane sheet. Similar results were re- 
cently reported in Refs@ andlsol. 

FigHUb) shows the Fermi surface evolution according 
to GGA-I-U calculations. In this case for AEp=+50 mcV 
the Fermi surface is similar to the GGA calculation, but 
shows a hole-like pocket near the (0,7r) point, whose ori- 
gin is the chain-apical band. This pocket merges with the 
Cu02 sheets at zero shift. This trend is again expected 
given the larger splitting of the chain-apical band at (0,7r) 
discussed in connection with Fig[T] For AEp=-55 meV, 
an electron-like pocket appears near (7r/2,7r), surrounded 
by a hole-like sheet. Going back to FiglU one immedi- 
ately realizes that this is also due to the enhanced split- 
ting in GGA-I-U: a similar pocket would appear in GGA 
for a much larger negative shift of over 200 meV. 

FigHKc) shows the PSIC results. The only structures 
in the Fermi surface are two large hole-like Cu02 sheets 
centered on (7r/2,7r). The Fermi level shift only moder- 
ately affect their area. No small pockets appear in this 
shift interval. 

Overall Fig[2] shows a marked sensitivity of the GGA 
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FIG. 2: Color on-line) Evolution of the Fermi surface of 
YBCO with Fermi-level shift AEp in the basal plane {k^ — 0). 
The main quantum oscillation orbits (-F„) are marked on the 
O meV and -55 meV panel for GGA (a), on the -1-50 meV 
and -55 meV for GGA-I-U (b). In PSIC there are no small 
pockets. 



and GGA+U Fermi surface to the relative positions of 
the bands. This suggests that subtle changes in dop- 
ing could result in the formation of small Fermi surface 
pockets. On the other hand, the PSIC Fermi surface is 
quite independent of doping, and would lead to predict 
or expect no small pockets at all. 



C. Fermi surface pockets: frequencies and masses 

To make contact with the quantum oscillations mea- 
sured in SdH and dHvA experimentSfii^i^i^ we calculate 
the quantum oscillation frequencies F={hA/2'Ke) from 
the cross-sectional area A of the orbits (i.e. the pockets), 
and the attendant cyclotron masses la^h"^ [dA/ dE) /2tt 
for the various structures found by the different tech- 
niques. They are reported in Figs. [21 HI and [5] for GGA, 
GGA-fU, and PSIC, respectively. 

For all techniques we report the high-frequency oscil- 
lations related to large cylinders; for GGA and GGA-I-U 
only, low-frequency oscillations related to small pockets 
are reported in a second panel. Thus, the frequencies 
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FIG. 3: GGA extremal dHvA frequencies (solid) for the large 
hole- like pockets (a) and for the small pockets (b). Dash- 
dotted lines are the calculated band masses. 
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FIG. 4: GGA-fU extremal dHvA frequencies (solid) for the 
large hole- like pockets (a) and for the small pockets (b). Dash- 
dotted lines are the calculated band masses. 



shown in Figl2{a), FiglUJa) and Fig\5\ (Fi and F2) are 
from the main Cu02 sheet surfaces, whereas those in 
FiglSJb) and Figd] (b) (F3 and F4) are from the smaU 
pockets. We note, first of all, that the frequencies calcu- 
lated for the main Cu02 sheets (Fi and F2) are similar 
for GGA and GGA-fU with frequencies between 3000 
T and 5500 T, whereas the frequencies calculated with 
the PSIC approach are between 1000 T and 2000 T. The 
reason of this difference is the lesser dispersion of the 
band structure as calculated with the PSIC technique. 
All values are way larger than the experimental one; the 
masses are typically a factor of two (or more) smaller 
than in experiment, and always negative. These Fermi 



surface sheets can therefore be ruled out as the origin of 
the experimental oscillations reported so far. 

Next we analyze the small-pocket signals in the fre- 
quency range to 900 T. In the GGA calculation, the 
hole-like pocket F3 has a frequency between 100 and 
600 T depending on the Ep shift; the experimental value 
would be attained at a shift of about -65 meV. The calcu- 
lated mass of this pocket is shift-independent, and equal 
to ~-1.4 me. The F4 pocket has a fairly low frequency 
of 100 to 300 T and a negative mass similar to F3. With 
the GGA-I-U approach we find the hole- like pocket F4, 
with a roughly shift-independent mass of ~-1.25 me and 
frequency in the 400 to 600 T range, and the electron-like 
pocket F3 with frequency between 400 and 800 T and a 
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FIG. 5: PSIC extremal dHvA frequencies (solid) for the large 
hole-like pockets. Dash-dotted lines are the calculated band 
masses. 



sharply varying mass, between 1.5 nig and 2.3 mg. 

Comparing with experiments, ^'^'^ several of our cal- 
culated pockets may seem good candidates. Frequencies 
and masses (in absolute value) are more or less in the ball- 
park. If we accept the assumption that the change of sign 
of the Hall resistance^ is purely due to the electron-like 
nature of the pockets, we implictly fix the experimental 
sign of the mass to a positive value. The frequency and 
mass deduced from observation^'^ would then be com- 
patible only with the F3 GGA+U pocket. 

A very recent measurement in YBa2Cu306.5i has re- 
vealed, in addition to the same signal of RefJ,, an oscil- 
lation with frequency and mass in the vicinity of 1600 T 
and 3.4 m.^ respectively. In Reffl (p-210) the signal is at- 
tributed tentatively to a hole- like pocket. In all our calcu- 
lations, including shifts, there is only one case (GGA+U 
at large negative shift, rightmost picture in panel b of Fig. 
[5]) in which hole and electron pockets coexist. Near the 
(7r/2,7r) point in the GGA-I-U calculations, starting at a 
shift of -55 meV, the structure recognizably involves two 
distinct pockets: one is the electron pocket F3 discussed 
above; the other a larger hole-like pocket surrounding F3 
itself. Their simultaneous presence is due to a change 
in curvature of the same band, most notably between 
(7r/2,7r) and (7r/2,0). The character of this band is, like 
that of F3, strongly chain-apical. The corresponding cal- 
culated frequency is about 2200 T and a mass of -1.4 
me. The frequency is very roughly similar to the 1600 
T measured in Ref|j, while the mass is over a factor two 
smaller. 

Overall, however, we conclude that there is not enough 



evidence to actually associate our calculated results to 
the experimental findings of Refs. il|l2t3|lJ. The reasons 
will be discussed in the next Section. 



D. Discussion and summary 

The calculations just reported have detected several 
small pockets (mainly hole-like) roughly compatible with 
the observed oscillation. However, all these small pock- 
ets have essentially chain or chain-apical character, and 
not in-plane character. GGA+U does seemingly finds 
the "right" pattern of coexisting electron and hole pock- 
ets, but (aside from the need for an artificial -60-meV 
Fermi level shift, corresponding to a 30% overdoping) 
both pockets have a chain-apical nature even stronger 
than the corresponding GGA-calculated band due to the 
remarkable (perhaps exaggerated) [/-induced lowering of 
in-plane Cu bonding states. 

On the other hand, there appears to be experimental 
evidence that the negative and oscillating Hall resistance 
at low temperature resulting from electron-like pockets 
(i.e. a positive mass) be related to states residing in the 
Cu02 planes. This is supported^ by the suppression of 
a6-plane conductivity anisotropy below 100 K, implying 
that chains do not conduct at low temperatures (and high 
field). 

Further supporting the fact that Fermi surface pockets 
are a plane-related feature, quantum oscillations were ob- 
served in YBa2Cu408.^'^'^ Calculations^iHiSS for that 
compound have shown that the GGA-calculated band 
related to the F3 hole pocket in YBa2Cu306.5 is now 
as far as 400 meV below Ep, hence cannot not rea- 
sonably invoked to explain the observations. Consis- 
tently, we found (unpublished calculations) that no pock- 
ets appear at all in the non-magnetic phase of chainless 

Yo.75Cao.25Ba2Cu306. 

We further recall that pockets appear only upon ap- 
preciably shifting the Fermi energy: the proper calcu- 
lated Fermi surfaces, i.e. those at zero shift, show no 
small pockets, except for the GGA F4 hole pocket of 
Fig. ^a) , related to the backfolding in the 2x1 cell of a 
pocket found by GGA itself in YBa2Cu307 (not seen by 
ARPES). 

Were we forced to embrace one of the methods applied 
here and the pertaining conclusions as the most reliable 
in this context, we would by all means pick PSIC, and 
conclude that in non-magnetic YBCO simply there are no 
small pockets, electron-like or otherwise. Indeed, among 
those used here, PSIC has shown to be by far the most 
dependable technique in the context of cuprates. For in- 
stance, the energy balance of various magnetic phases 
of YBa2Cu306+a; is correctly described, and so are the 
general properties of a number of cupratesJ^ii^ii^ Fur- 
thermore, in the context of Fermi surface determination, 
PSIC matches ARPES perfectly for YBa2Cu307 (unpub- 
lished calculations) whereas GGA finds, as mentioned, a 
zone-corner pocket which ARPES does not observe. 
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In summary, we presented calculations of the elec- 
tronic structure of YBCO in the non-magnetic state with 
three different DFT-based approaches: GGA, GGA-I-U 
and PSIC. Upon substantial shifts of the Fermi energy, 
GGA and GGA-I-U do produce small Fermi surface pock- 
ets, mostly originating from chain or chain-apical bands, 
with frequencies and band masses similar to those exper- 
imentally observed (one GGA-f U pocket has a positive 
cyclotron mass, i.e. is electron-like), while PSIC shows 
no small pocket at all. As discussed, our conclusion is 
that there is no unambiguous evidence for the existence 
of electron-like pockets -nor, indeed, of any pockets- in 
the non-magnetic metallic state of YBa2Cu306.5. In ad- 
dition, no pockets (either electron or hole) derive from 
in-plane states. This is a conclusion coherently obtained 
by three different ab initio techniques. We suggest that 



the experimentally observed pockets are a property of 
another state of YBCO, possibly characterized by some 
form of ordering (probably magnetic, given its coex- 
istence with superconductivity up to high doping re- 
vealed by many experiments) causing a Fermi surface 
reconstruction. We will present elsewhere further first- 
principles work in this direction. 
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